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The calculation of nonnormal tropospheric bending of radio waves is treated in terms 
of a reduced-to-sea-level value of the refractive index. This method emphasizes departures 
of bending from the average bending for the United States and consists of visualizing ray 
bending as consisting of two parts; an "average" component and a "departure-from-average" 
component. The "average" component comprises most of the bending and is obtained 
accurately from refraction tabulations while the component due to departures is easily 
obtained by graphical means. 



1. Introduction and Background 

The angular bending, r, of a radio ray passing from 
the earth's surface where the refractive index is 
n s to any point in the atmosphere of refractive 
index n is given by 
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where 6 is the local elevation angle of the ray. (The 
geometry is shown on fig. 1.) In practice, (1) is 
evaluated by numerical integration, assuming that 
n is stratified in layers concentric with the earth and 
decreases linearly between measured or assigned 
values of n [1, 2]. 1 Although at any instant the 
vertical distribution of n may be quite complex, it 
has been customary [3] to assume that n decreases on 
the average linearly with height up to ){ km, thus 
vastly simplifying many practical problems en- 
countered in radio engineering applications involving 
antennas at these heights. Even though the assump- 
tion of a linear decrease of n is adequate for the 
majority of applications, the advent of space com- 
munications has forced a closer consideration of 
actual atmospheric n structure [4, 5, 6] with the 
general conclusion that it is more realistically repre- 
sented by an exponential decrease with height than 
by a linear decrease. Millington [7] and Norton 
[8, 9, 10] have given numerous examples of the use 
of these more realistic atmospheres in the calculation 
of transmission loss in both ionospheric and tropo- 
spheric propagation, while Wait [11] has given rigor- 
ous series formulas for the distance to the horizon 
in such an atmosphere. 

In fact, the International Radio Consultative 
Committee has recently recommended for inter- 
national use a basic reference atmosphere to describe 
the average refractive index structure of the at- 
mosphere by 



* Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Geometry of atmospheric refraction of radio waves. 



N(h) = [n(h)-l] 10 6 =iV s exp(-ch). 



(2) 



Tables of refraction variables have been prepared 
[12] for atmospheres of this form for specific values 
of N 3 and c throughout their normal range of 
variation. 

The systematic exponential decrease of the re- 
fractive index with height is relatively large compared 
to variations introduced by fronts and airmass dis- 
continuities. To obtain a clear picture of such 
variations one should eliminate the systematic 
height changes from the actual observed value, 
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N(h), and thus emphasize the departures from normal 
of the actual N structure. This may be accom- 
plished by 



A(h, N s )=N(h)+N s [l-exp(-ch)) 



(3) 



which has been found [13] to be more effective than 
previous corrections involving the assumption of a 
linear N(h) distribution [14 through 17]. 

The purpose of this paper is to consider departures 
of n structure from normal, as specified by (2), the 
effect of these departures upon ray refraction, and 
to present a new method for calculating r based 
upon these departures which will differ from other 
approximate methods already in the literature [1, 
2, 5, 12, 18, 19, 20]. 

2. Calculation of Nonnormal Refraction 

In evaluating the bending of a ray one is concerned 
with the refractive index gradient: 



dA(h, N s )=dN(h)+N s c exp(-ch)dh 



(4) 



or 



10 6 dn=dN(h)=dA(h, N s )-N s c exp{-ch)dh. (5) 

The integral for the bending of radio waves, 
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The height variation of 6, although dependent upon 
the refractive index distribution through Snell's 
law, is at least four times more sensitive to purely 
geometric changes than to refractive index changes 
[2]. This then suggests that we may simplify our 
problem by assuming that 6 may always be deter- 
mined by its value in the N(h)=N s exp(-ch) at- 
mosphere for the height under consideration. When 
this assumption is made the second term on the 
right-h and s ide of (6) becomes simply the bending 
in the N(h)=N s exp(-ch) atmosphere, designated 
r(h, N s ). Equation (6) may now be written 



r l,2= — 
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dA(h,N s ) + r(h,N s ). (7) 



Since r (A, N s ) normally comprises most of the bending 
and is accurately tabulated (for several specific 
values of N s ), approximate or graphical methods 
may be used to determine 
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The perturbation term represented by eq (8) pre- 
sumably contributes only a small portion of r which 
is taken as justification for the further approxima- 
tions that n=\ and that cot 0=1/0, then 
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The problem now consists of determining the best 
method of estimating r(h, N s ) and the variation of 
with height. One method would be to refer the 
method to a single exponential N profile typical of 
the United States : 



iV(A)=313exp (-0.144/*,) 



(10) 



where h is in kilometers. Equation (10) represents 
average conditions and would have the advantage 
of referring all A 7 structure to a single "standard" 
atmosphere. However, climatic and synoptic 
changes in refractive index profiles can be accounted 
for by utilizing the range of refraction variables 
tabulated in reference [12]. These tabulations are 
in terms of specific values of N s ranging from ^=200 
(extremely dry and high locations) to ^ = 450 
(extremely hot and humid locations). However, 
refraction effects are more sensitive to initial N 
gradient conditions than to the value of N $ [5], 
particularly at small elevation angles. Thus, in 
what follows, the choice of a reference atmosphere 
will be determined by the initial gradient of the N 
profile for which we wish to determine the bending. 
Table 1 lists the exponential profiles for which 
refraction tabulations are available. The N s value 

Table 1. Initial N gradient J} AN E , in the C.R.P.L. exponential 
reference atmosphere ° 



Range of A N t (2V units/km) 



- AN, < 27. 55 
27.55<-AiV e <35.33 
35.33<-AA T e <42. 13 
42.13<-AiV e <49.52 

49.52<-AN 6 <59.68 
59.68<-AiV £ <71.10 
71. 10<-AiV e <88. 65 

88.65<-AiV € 



N(h) 



200 exp(-0. 1184/i) 
252.9exp( -. 1262/0 
289 exp( -.1357/i) 
313 exp( -. 1438A) 

344.5 exp( -.1568ft) 
377.2 exp( -. 1732/0 
404.9 exp( -.1898/1) 
450 expf -.2232ft) 



a Note: height, ft, is in kilometers. 

and the range of initial N gradients assigned to 
each atmosphere are also listed. Since, theoretically, 
there are an infinite number of exponential atmos- 
pheres but, in practice, onh T a few with tabulated 
refraction variables, it will be necessary to assign 
each actual JV profile encountered to a specific 
exponential reference atmosphere according to its 
initial N gradient. The initial N gradient of 



N(h) = N s exp (-ch) 



is given by 
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\ dh Jo 



(11) 



For any actual case, however, one would not be 
able to define (dN/dh) since the radiosonde reports 
at only certain points above the ground and not 
continuously. Thus if the first reporting level is 
at a height of e, the actual initial gradient is deter- 
mined from 



-AAT £ : 



N.-N. 



(12) 



For the practical application that is the objective 
of this paper we shall approximate (dN/dh) by 
AN e , and the choice of an exponential atmosphere 
from table 1 will be based upon the observed value 
of AN e . As a consequence of this choice N s will 
usually differ from A T (0) . 

The notation A(h, 313) and r(h, 313) indicates, 
for example, that the initial gradient of the particular 
profile under study falls into the 313 exponential 
atmosphere category, and that A, r, and 6 are 
determined from that particular atmosphere. Note 
that .4(0, 313) is not equal to 313 but A(h, 313) will 
approach 313 at sufficiently large heights. 

Radiosonde data from Santa Maria, California 
(a good test station because the refractive index 
profile departs quite radically from any smooth 
model) were graphically converted into A(h, 313) 
units (see fig. 2) and the bendings calculated for 
several values of O . The A(h, 313) profile is com- 
pared on figure 3 with the departure of the bendings 
from the values expected in the 313 exponential 
atmosphere and emphasizes that r is an integrated 
effect of all refractive index departures from the 
surface to the point under consideration while 
A(h, 313) is determined solely from data at the 
height h alone. For instance, a strong surface- 
based refracting layer will bias all bending values 
above the layer while A(h, 313) is independent of 
the values of N at any other height. It is quite 
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Figure 2. N 8 [1-exp (-ch)] versus height for conversion of 
N(h) profiles into A(h, N 8 ) profiles. 
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Figure 3. Comparison of A(h, 313) with the departures of 
ray bending from the value in the 313 exponential atmosphere 
for Santa Maria, Calif., Sept. 4, 1952, 1900 local time. 

Note that negative gradients of A(h, 313) produce positive values of ray bending. 

clear that the agreement in height structure between 
the bending departures and A(h, 313) is very pro- 
nounced on figure 3, especially in the layers nearest 
the ground and for small values of O , indicating 
that a knowledge of A(h, 313) alone could enable 
one to obtain a visual impression of the effect of 
nonnormal atmospheric structure upon the refraction 
of radio waves. 

The above illustration has served to introduce a 
method of calculating r in terms of departures of A^ 
structure from normal. The following section will 
be concerned with testing the accuracy of the 
method. 

3. Comparison of Bendings 

Significant level data for the radiosonde observa- 
tions taken at 11 geographically and climatically 
diverse U.S. Weather Bureau observatories were 
analyzed to give adequate coverage of climatic con- 
ditions found within the com try. Individual radio- 
sonde observations taken during the years 1951 to 
1952 were examined for each observing site, and 
profiles were selected to give a measure of the maxi- 
mum variation of bending due to both geographic 
and local refractive index profile variations. The 
bending was then calculated for each of these profiles 
by both the standard numerical integration of (1), 
[1], and the departures-from-normal method de- 
scribed above. The differences between the bend- 
ings obtained by the two methods are given in table 

2 as a percentage of the value obtained by conven- 
tional methods. A classification of the data of 
table 2 was made as regards height. Since most of 
the total bending to 70 km occurred within the first 

3 km [5], the integral correction involving A was 
only applied to that height, the bending in the 313 
atmosphere being used to estimate the average 
bending above 3 km thus assuming that the effect 
of iV departures above this height is negligible. The 
mean errors are quite close to zero, indicating the 
absence of any marked bias in the method of calcula- 
tion. One also notes from table 2 that the maximum 
percentage difference occurs at 00=0 and generally 
decreases with increasing initial elevation angle. 
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Table 2. Percentage deviation of ray bending values calcu- 
lated by the departures-from-normal method from convention- 
ally calculated values 



0o(mr) 



ft (km) 



Mean % of error 

Std deviation of % errors. 
Max % error 



Mean'%*of error 

Std deviation of % errors. 
Max % error 



10 



52.3(3°) 



A. Bendings calculated relative to multiple 
exponential atmospheres 



-0.84 

3.19 

-10.0 


-0.35 
3.65 

+8.6 


+0.49 
3.59 

+9.4 


+0.30 

3.71 

+9.4 


-0. 007 

.97 
+2.5 


-0.80 
3.95 

+8.5 


B. Bendings calculated relative to the 313 exponen- 
tial atmosphere 



-0.03 


+0.04 


-3.04 


-1.14 


+0.78 


7.57 


7.15 


6.36 


5.19 


1.14 


14.5 


-14.0 


-14.6 


+10.9 


-2.7 



-0.87 

4.29 

+9.0 



Since the mean deviations are less than 1.0 percent 
the present system of calculation appears to be 
adequate for normal applications. One would prob- 
ably not wish to use it for estimating t for O ^>52.3 
milliradians (mr) since at such high angles the total 
bending is as accurately estimated solely from 
ground level values by [18]: 



r=iV s X10- 6 cot d . 



(13) 



One might also wish to use the departures-from- 
normal method referenced to a single average N 
profile t3 r pical of the entire United States. This 
approach was investigated by use of such an average, 
the 313 exponential atmosphere, and the results are 
tabulated in table 2. It is quite evident that the 
use of multiple reference profiles, classified according 
to AN e , produces significantly smaller standard 
deviations and more consistent mean errors than 
does the use of a single average profile. This is, of 
course, what one would expect since the AiV € method 
adjusts the profile to the initial la}^er where a signifi- 
cant percentage of the total bending occurs. 
Interestingly enough, the maximum errors for 6 =0 
arise from profiles such as shown on figure 2 where 
a very strong elevated layer exists well above the 
surface layer. Under this condition, the classifica- 
tion of profiles according to initial gradient does not 
account for the major anomaly of the profile. The 
maximum errors for 6 = 10 mr are due to a quite dif- 
ferent cause; the presence of ducting profiles in the 
test profiles. Under ducting conditions rays may 
be traced only for initial elevation angles greater 
than the angle of penetration [21] which, for radio- 
sonde observed ducts, is less than 6 mr [22]. For 
this reason the data for O =1O mr includes many 
profiles that represent near ducting conditions and 
consequently maximum departures from normal re- 
fraction. In fact, the maximum errors for O = 10 mr 
are associated with profiles having very thin surface 
ducts with normal or subnormal profiles above the 
duct, a condition that produces overestimation of 
the bending above the duct. 



It should be observed that the A(h, N s ) units as 
used in this paper are always referenced to zero 
height above ground level. We are able to ignore 
the altitude of the ground level by using multiple 
reference atmospheres which have, in effect, a 
built-in dependence of N s upon ground level alti- 
tude. This effect may be accounted for with a 
single exponential atmosphere by adding the alti- 
tude of the ground level, b, to the height, h. Thus 
the equation for bending becomes 



-r 



n * 10~ 6 cot 9 



n 



dA(h, 313) + exp (-c&)r (A, 313). 

(14) 



The exponential correction factor to r(h, 313) is a 
constant for any particular station and is important 
only when b is large and consequently N s appreciablv 
less than 313.0. 

One should bear in mind, however, that the above 
errors are about one-half that of other simple methods 
of estimating the bending of radio rays for 0o<C52.3 
mr such as regression techniques [18] or model 
atmospheres based upon the surface value of initial 
gradient conditions [5]. It would seem then that 
an additional advantage of the present system lies 
in the facility witli which the user would come to 
appreciate the relative importance of various profile 
departures. 

The data required for use of the present method 
are contained in reference [12]. Values of 
N s [l-exp(-cA)] are given in table 3. 

Table 3. Values of N 8 [1-exp (-ch)] 



Ht. 


200.0 


252.9 


289.0 


313.0 


344.5 


377.2 


404.9 


450.0 


(km) 


















0.01 


0.2 


0.3 


0.4 


0.4 


0.5 


0.6 


0.8 


1.0 


.02 


.5 


.6 


.8 


.9 


1.1 


1.3 


1.6 


2.0 


.05 


1.2 


1.6 


2.0 


2.2 


2.7 


3.2 


3.9 


5.0 


.10 


2.4 


3.1 


3.9 


4.5 


5.3 


6.5 


7.6 


9.9 


.20 


4.7 


6.3 


7.7 


8.9 


10.6 


12.8 


15.1 


19.7 


.50 


11.5 


15.4 


19.0 


21.7 


26.0 


31.3 


36.6 


47.5 


1.00 


22.3 


30.0 


36.7 


41.9 


50.0 


60.0 


69.9 


90.0 


2.00 


42.2 


56.4 


68.7 


78.3 


92.7 


110.4 


127.6 


162.1 


5.00 


89.4 


118.4 


142.4 


160.5 


187.2 


218.6 


247.8 


302.6 



4. Conclusions 

The present study has shown that deviations from 
normal ray bending are mirrored by departures of 
refractive index from values expected in average 
atmospheres of exponential form. Further, this 
departures-from-average method emphasizes that 
^-profile anomalies cause the greatest change in the 
bending when they occur near the ground, these 
changes being larger the smaller the initial elevation 
angle. A graphical method based upon these 
departures may be used to calculate the bending of 
radio rays. 

The authors extend their gratitude to K. A. 
Norton and G. D. Thayer for their many critical 
and stimulating reviews of this work. 
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